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Abstract 

Based on Keldysh's nonequilibrium Green function method, the spin-dependent transport prop- 
erties in a ferromagnet-quantum dot (QD)-ferromagnet coupled system are investigated. It is shown 
the spin current shows quite different characteristics from its electrical counterpart, and by chang- 
ing the relative orientation of both magnetizations, it can change its magnitude even sign. The 
current-induced spin transfer torque (CISTT) is uncovered to be greatly enhanced when the bias 
voltage meets with the discrete levels of the QD at resonant positions. The relationship between 
the CISTT, the electrical current and the spin current is also addressed. 

PACS numbers: 75.47.-m, 73.63.Kv, 75.70.Cn 
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I. INTRODUCTION 



Investigations on the spin-dependent transport in a magnetic tunnel junction (MTJ) 
have attracted much interest in the last decade, as the MTJs have essential applications in 



spintronic devices jl, 0, 0, U, 0]. One of phenomena in MTJs is the so-called tunnel mag- 
netoresistance (TMR) effect, which states that the tunneling current through the junction 
depends sensitively on the relative orientation of the magnetizations of both ferromagnetic 
electrodes, that is caused by the spin- dependent scattering of conduction electrons. The dif- 
ference between the currents through spin up and down channels in MTJs is usually referred 
to as the spin current yi. On the other hand, a reverse effect to TMR was predicted indepen- 
dently by Slonczewskip] and BergerjTj], i.e., the spin-polarized electrons passing from the left 
ferromagnetic layer into the right layer in which the magnetization deviates the left by an 
angle may exert a torque to the right ferromagnet and can change the orientation of its mag- 
netization. This effect was coined as the spin transfer effect (for review, see e.g. Ref . jl], 0] ) , 
which can lead to the current-induced magnetization reversal, and might offer a promise for 
the current-controlled spintronic devices. It has recently attracted intensive investi gati ons 



both experimentally 
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Among other magnetic tunnel structures (e.g. Refs.[23j), the ferromagnet-quantum dot- 
ferromagnet (FM-QD-FM) coupled systems have also received much attention recently. Pre- 
vious theoretical works on spin-dependent transport through QDs are ma inly focused on the 



tunnel electrical current and the TMR effect for co 
and noncollinear configuration 
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linear configuration 
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391 ]. However, the inves- 



tigations on the spin current and the current-induced spin transfer torque (CISTT) in such 
systems are sparsely reported. It is the purpose of this paper to study the spin current and 
spin transfer torque in FM-QD-FM coupled MTJs. 

In terms of Keldysh's nonequilibrium Green function method, the tunnel electrical cur- 
rent, the spin current and the spin transfer torque in the FM-QD-FM coupled system will 
be investigated. It is found the spin current exhibits different behaviors from its electrical 
counterpart. It is also shown that the resonant positions for the tunneling electrical current 
and spin current become far separated with the increase of the Coulomb interaction U in 
the QD. The magnitudes of the CISTT at two resonant positions are found to be greatly 
enhanced at resonant positions. At the bias larger than (eq + U)/e, the CISTT can reach 
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a saturation plateau which is independent of the Coulomb interaction. The CISTT shows 
a kink-like behavior with the increase of the spin current. The relationship between the 
CISTT, the electrical current and the spin current is also addressed. 

The rest of this paper is outlined as follows. In Sec. II, the model and formalism for the 
tunnel electrical current, the spin current and the current-induced spin transfer torque will 
be established. In Sec. Ill, the corresponding numerical results will be given. Finally, a 
summary will be presented. 

II. MODEL AND FORMALISM 

Let us consider a single-level QD coupled to two ferromagnetic electrodes. The left (L) 
and right (R) electrodes are connected with the bias voltage V/2 and —V/2, respectively, as 
shown in Fig. 1. The magnetic moment of the left FM is assumed to be parallel to the 
z axis, while the moment Mr of the right FM is aligned along the z' axis which deviates the 
z axis by a relative angle 9. The tunnel current flows along the x axis and perpendicular to 
the junction plane. The system consists of the left, right, QD and coupling parts, and can 
be described by the following Hamiltonian: 

H = H L + H R + H d + H T , (1) 

with 

Hh = ^2 £kLaa kLa a kLv, (2) 
ka 

H R = y^Jenjk) - aM R cos 9]al Ra a k Ra ~ M R sin 9a[ R(T a kR ^, (3) 

ka 

H d = y^epclcg + Un^ni, (4) 

a 

H T = y^Tkaa\ aa c a + h.c, a = L,R. (5) 

kaa 

where e kaa = e a (k) — <rM a — eV a is the single-electron energy for the wavevector k and 
spin a in the a electrode, a kaa and c a are annihilation operators of electrons with spin a 
in the a electrode and the QD, respectively, n a = cJ.c CT , U represents the on-site Coulomb 
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FIG. 1: Schematic illustration of the system consisting of two ferromagnets and a quantum dot 
separated by the tunnel barriers, where T^ a (a = L, R) stands for the coupling matrix between the 
a electrode and the QD, and both magnetizations are aligned by a relative angle 6. 

interaction between electrons in the QD, and T^ a is the coupling matrix elements between 
the a electrode and the QD. 



A. Tunnel electrical current and spin current 

The tunnel electrical current is composed of the sum of the currents carried by spin up 
and down electrons: 

I(V) = I L1 (V)+I Ll (V), (6) 

while the spin current is defined by the difference between the electrical currents through 
the spin up and down channels jj]: 

I s (V) = I L ^(V)-I Ll (V), (7) 



where 



h 

kL 



h 

kL 



with G^^it.t 1 ) = i(a\. La (t')c a '(t)) the lesser Green function. 
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By applying the Langrenth theorem and Fourier transform, we may obtain the following 
equation 

GUe) = n L [G r (e)g< L (e) + G<(e)g a kL (e)}, (8) 
where G^s) = [ V G^s) = ( ^ ^ | , ,<(*) 



i27r/ a (efe a t)5(e - £tet) 

i27rf a (e ka i)S(e - e^) 



i 



i 



e-e kal ±ir] 

G r (e) and G < (e) are the retarded and lesser Green functions of the QD, respectively, and 
f a (e) is the Fermi distribution function in the a electrode. By defining M(e)= [fL(£)(G r (e) — 
G a (e)) + G < (£)]T L (e), we can simplify Ift$(V) as 

zc f els 

I m (V) = -- I -M mi) (e), 

where M = ( | is a 2 x 2 matrix, T a (e) = \ \ with T aa (e) 

2nY J \T ka \ 2 S(e - s kaa ). 

ka 

The lesser Green function G K can be calculated by the Keldysh equation G < = G 7 "E < G a . 
To get E<, we invoke Ng's ansatzQ: S< = S<5, where £<(e) = i[f L (e)T L + f R (e)TZT R TZ% 
B = (TT - Eg)" 1 ^ - E5(e) - Eg(e) = -i[r £ + ^r^t], 57(e) - = G" 1 " 1 - G""~\ 
cos | — sin | 

and 7Z = \ " . Under these considerations, we finally get 




sin | cos ~ 



%)^) = |/ ^(/«"/^TT(U)- (9) 



where X = G r KT R lV BG a T L = | u | is also a 2 x 2 matrix. 



Consequently, the tunnel current and spin current have the forms of 

I(V) = C -j ^(f R -f L )TrX, (10) 
UV) = e -f ^-(f R -f L )Tr(Xa 3 ). (11) 

,10. 

with £T 3 = I | the Pauli matrix. 

-1 
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B. Current-induced spin transfer torque 



The spin torque exerting on the right ferromagnet is determined by the time evolution 
rate of the total spin in the right ferromagnet 0, [J , which is composed of two parts: one is 
caused by the spin-dependent potential that is known as the equilibrium torque, the other is 
the current-induced spin transfer torque caused by the tunnel Hamiltonian Ht- By means 
of the none quil ibrium Green functions, the CISTT exerting on the right ferromagnet can be 
obtained by|l7j 



Br' 

T ux = - C os 



me H [ ^[^(e)^] + sin^e^T f ^-Tr[G< R (e)a 3 T* R ], (12) 



kR " kR 



where &i 



1 

1 

By a treatment similar to Eq. (jSji . we can obtain 



GUe) = T kR [G r (e)Kg< R (e)ri + G<(e)TZg a kR (e)^}. 

Therefore, the CISTT can be rewritten as 

t Rx ' = ^ J de{f R -f L )Tr[G\e)Y L BG a {e)nY R n\-cos6a 1 + smea 3 )\. (13) 

The remaining task is to calculate the retarded Green function G r . 
By the equation of motion, we can derive 

(e - e )G r (e) = 1 + Y?LG r ka (e) + UG^(e), (14) 

ka 

( G^ r (e) G n ^ r (e) \ 

where 1 is a unit matrix, G^ r (e) = , GZ^ r (t - f) = -iO(t - 

\GM>(e) G^>(e) J 

f)({aw(t),4(f)}>, G°'°mr(t- t >) = _ id (t - t ')({cAt)n*,(t),cUt')}), and 

g \l = T kL g r kL G r '(e) , 
G kR = T kR ngl R rtG r (e). 

Substituting these equations into Eq. ([14)1 . up to the third-order of Green functions, we 
arrive at 

(e - e - U)G^'(e) = N + ^[T^G^e) + T ka G 2 k ^(e) - T^G 3 ^ \e)\, (15) 

ka 



\ n l) -\ C | C T/ I n i(2)r _ I ^ka U kc 



G a k ':' mr (t-t') = -i9(t-t')({a ka At)nAt),cUt')}), 

Gt:' mr (t-f) = -•e(t-o<{oL'(*)^(*)^(*).4(0}>, 

G ^,mr {t _ f>) = _^( t - t ')({aw(t)4,(t) C(T ,(t), C t(t')}). 



k'a 
k'a 

dl^Tis) = T kL [G r (s) - G^(e)) + £[-T* a C$£ a ( e ) + T* a Gt^ a (e)], 

k'a 



with ~f ka (e) = ( — , I = I — — ^ x | , 



by: 




i 



' ^ka * 

£— 2eo— C+efc Q iT+* ? / 

and the third-order Green functions defined 



r ,cr'o-,l(3)r /, 
^kLk'a \ L 





= -i9(t 


-t'){{a kLa >{t)a\, a ^{t)c ff 


it), 


cl(t')}}, 


(~<cr'cr,2(3)r /. 
^kLk'a \ l 





= -i6(t 


- t'){{a kL(Tl {t)cy,{t)a k t a& 


(t), 


cl(t')}}, 


^cr'a,3(3)r /, 
u kLk'a \ l 





= -%6(t 


-t')({al L& ,(t)a k , aa r(t)cv, 


(t), 


cum, 


r a'crA(3)r ( 
^kLk'a \ L 





= -i9(t 




it), 


4(01), 


(~<cr'cr,5(3)r /. 
^kLk'a \ l 





= -i9(t 


-t'){{akLa'{t)a\, aS ,{t)c & 


it), 


cl(t')}}, 


x-ycr'(T,6(3)r / , 
U kLk'a \ l 





= -i6(t 


- t')({akLa>{t)cl,{t)a k > a& 


it),ct(t')}). 



It is worthily mentioning that one usually takes the Hatree-Fock decoupling approximation 
up to the second-order for Green functions in most systems. Such a decoupling scheme 
is not adequate for the present single-level QD coupled system, as it could smear some 
characteristic features that come from electronic correlations. In order to extract more in- 
formation from the many-body interactions, we should consider the equation of motion up 
to the third-order of Green functions G l j^ r (e) in Eq. ()15p. In above derivations, we have 
invoked the following decoupling approximations for Green functionsj^: G^y^ r (e) = 
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G kLWa )r ( £ ) = 8kL,k'a[l ~ j liskLa'^G' 7 '(e) .As we ignore the spin-flip scatterings in the 
present system, the off-diagonal elements associated with different spins of the third-order 
Green function are sent to zero, namely, the higher-order spin correlations in the ferromag- 
netic leads are neglected. 

With these results, we may obtain 

GT(e) = T kL9kL G^(s), 

G 2 k { f\e) = n L f ka [G^(e)-F L G r (e)}, 

GfT{e) = -T kL gl a [G^(e)-F L G^e)}, 



fa(e kal ) 



with F n = | I and F a 

Jackal) J y /a(£feat) 

Similarly, a more straightforward but somewhat complicated calculation gives rise to the 
following equations 

1, 
T 

+T kR 



GlT(e) 



Gfi r (e) 



2 ^ r / \ • 2 ^ r / \ 

cos -g kR (e)+sm -g kR {e) 



(16) 



r 2(2)r 
^kR 



-T kR ^n9[giUe)f R (e km ) - Slfte)/*^)]*!^) 



-T kR cos 2 -g r kR (e)F R + sin 2 -g r kR (e)F R 



G r {e) 



+T kR 



cos 2 d -9 r kR {e) + sin 2 ^g r kR {e) 



G^ r (e), 



(17) 



where U\ 



GlT(e) 



-T, 



kR 



kR 



2 ^ r / \ 2 ^ r / \ 

sm ^9 kR {e) + cos -g kR (e) 



9 6 
sin 2 ^9 r kR (e)F R + cos 2 -g r kR (e)F R 



G (2)r {e) 

G r (e), 




and gl a (e) 



~gl r a (e) 



0fca( e ) 



(18) 



e-2eo-U+e ka ^+iri 



1 



e-2s -U+e ka i+iri 
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On the other hand, the second-order retarded Green's function can be expressed as of 
the following form 

W 1 G (2)r (e) = N + W 2 G r {e), (19) 



where 



kL 



kL 



kL 



.2 ® r 



W x = (e - £ - U) - YTkLT* kL g r kL {e) - ^TkT^^) - jy kL T^ L gl L (e) 

~ Yl TkRT k 

kR 

-YT kR T* kR ^ sm9 tf£(e) - §&(e) 

kR 

J-kR-l-kh 



0\ 



kR 



cos 2 d -g r kR {e) + sin 2 9 -gl R {e) 



1 kRJ- k R 



kR 



cos 2 ^g r kR (e) + sin 2 ^g r kR (e) 



W 2 = -YTkLTt L ~gl L {e)F L - ^T^^e)^ 

kL kL 

+Y^T kR T* R ~ sin e[gl^(e)f R (e kR1 ) - gi^Ue^a, 
kR 

+X;rjw%iaine[^(e)/H(e^)-^(e)/ il (e J ^)]^ 

kR 

Erri rri* 
1 kR-L kL 



kR 



9 9 
cos 2 -g r kR (e)F R + sin 2 -g r kR (e)F R 



1 kRJ- k R 



kR 



9 9 
cos 2 ^g r kR (e)F R + sin 2 -g r kR (e)F R 



By combining Eq. (|19p with Eq. f|15j) . we get 



(e - e - S r - UW{ l W 2 )G r = 1 + [Ay^iV. (20) 

From this equation, the retarded Green function G r can be obtained. We would like to point 
out that the averaged values involved in Eq. (fT5j) should be obtained self-consistently by 



de 
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FIG. 2: The bias dependence of the differential conductance G for different polarization P in 
parallel (a) and antiparallel (b) configurations of magnetizations. The parameters are taken as 

e = ior , U = 40r , k B T = o.3r . 

These above-mentioned equations establish the fundamental basis for numerically investi- 
gating the spin-dependent transport properties of the FM-QD-FM coupled MTJs. 

III. NUMERICAL RESULTS 

Without losing the generality, in the following numerical calculations we may further 
suppose that the two ferromagnets are made of the same materials, i.e., Pl = Pr = P, 
where P L{R) = (T Lm - r L(i?)i )/(r L(i?)T + T L{R)l ) is the polarization of the left (right) 
ferromagnet. Then, we can introduce IV^ = r^f^ = r (l ± P), where r = T L ( R ^(P = 
0) = Tl(r)i(P = 0) will be taken as an energy scale. The tunnel matrix elements Tki(R) is 
presumed to take values at the Fermi level. We will take Io = ^ and Go = ^ as scales for 
the currents and corresponding differential conductance, respectively. 
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A. Differential conductance 



The bias dependences of the tunnel electrical current and the differential conductance 
G = in the FM-QD-FM coupled system have been reexamined within the present scheme. 
It is found that the tunnel electrical current exhibits two-step features, which corresponds 
to the resonant tunneling of electrons through the QD at energy levels Eq and Eq + U. The 
corresponding differential conductance as a function of the bias thus exhibits resonant peaks, 
where the large peak at a lower voltage corresponds to the discrete level Eq, and a small 
resonant peak corresponds to the level Eq + U. It is uncovered that with the increase of the 
interaction U, the peaks of G(V) corresponding to the discrete level Eq + U shift towards the 
higher voltage side without changing the peak amplitude for both parallel and antiparallel 
configurations. This is because an increase of U makes two discrete levels Eq and Eq + U 
become far apart, thereby leading to that the resonant positions of main and charging peaks 
of G(V) are far separated. These observations are quite consistent with the previous studies 
(see e.g. Ref.Q). 

The bias dependence of the differential conductance for different polarizations in the 
parallel and antiparallel configurations of magnetization is shown in Fig. 2. It can be found 
that in the parallel configuration there is only slight changes for different polarizations except 
for P = 1, as displayed in Fig. 2(a). We can understand this property from the point of view 
of the resonant tunneling through two barriers|31|. namely, the two barriers are of the same 
height, and the resonant probability is unity in this case, that is independent with P. The 
small changes with P observed in Fig. 2(a) is caused by the interaction U. However, it is 
quite different for the case with full polarization P = 1, which comes from the fact that the 
conductance for P = 1 has only one channel. It is interesting to note that with increasing 
P in the case of parallel alignment, although the main peaks at the bias corresponding to 
the energy level Eq are suppressed, the charging peaks at the bias corresponding to the level 
eq + U are enhanced. In this case, more electrons from the left enter into the QD through 
the channel Eq + U to tunnel into the right lead. For the antiparallel configuration, it can be 
seen that the main peaks of the conductance decrease as the polarization P increases, while 
the charging peaks remain almost intact, as shown in Fig. 2(b), that is different from the 
case of parallel configuration, and is due to the conventional spin- valve effect in MTJs. For 
P — 1, the conductance becomes zero, which is nothing but a perfect spin- valve effect. 
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FIG. 3: The bias dependence of the spin current I s (a) and spin differential conductance G s (b) 
for different angle 9 at P = 0.4. The other parameters are taken the same as in Fig. 2. 

B. Spin current and spin differential conductance 

The bias dependence of the spin current I s and its corresponding spin differential con- 
ductance defined by G s = ^ for different angles is shown in Fig. 3. It can be seen that, by 
changing the relative angle 9, the spin current can change its magnitude even its sign. At 
V > 0, when 9 < n/2, the spin current I s shows step-like behaviors, and I s > 0; while for 
9 > 7r/2, I s exhibits behaviors similar to a basin-like shape, and I s < at the bottom of the 
"basin". This shows that the effect of the relative orientation of magnetizations on the spin 
current is more obvious. At two resonant positions corresponding to levels £q and Eq + U, I s 
begins to sharply change. It is clearly seen in Fig. 3(b) that at V > 0, the spin conductance 
shows sharp resonant peaks and dips at the levels e and e + U, although for 9 = and 
7r/3 the resonance corresponding to Eq shows a peak, while for 9 = 7r/2, 27r/3 and 7r the first 
resonance shows a sharp dip. The resonances corresponding to the level Eq + U exhibit peaks 
for all 9. At 9 = 0, the two resonant peaks of I s have comparable weights; at 9 = 7r/3, the 
first peak has a small weight, while the second peak has a large weight; at 9 = tt/2, the first 
dip of I s has a small amplitude, while the other peak has a large amplitude; at 9 = 27r/3, 
the resonant dip and peak have almost the same amplitude; and at 9 = n, the first dip has 
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FIG. 4: The bias dependence of the spin current I s for different polarization P in parallel (a) and 
antiparallel (b) configurations of magnetizations. The other parameters are taken the same as in 
Fig. 2. 

a larger amplitude than the peak. Evidently, the spin current and spin conductance have no 
spin- valve effect, showing that the spin current and spin conductance have quite different 
characteristics from their electrical counterparts. 

The step-like and basin-like behaviors of the bias dependence of the spin current for 
different polarizations in parallel and antiparallel configurations are presented in Figs. 4(a) 
and (b), respectively. For P = 0, as there is no polarization in electrodes, the contribution of 
electrons with spin up is the same as that of spin down, leading to the spin current vanishes. 
In the case of parallel alignment (8 = 0), with increasing the polarization P, the contribution 
of spin up is increased while that of spin down is decreased, thus resulting in the magnitude 
of the spin current increased, as shown in Fig. 4(a). The step-like behavior of I s against V 
comes from the resonant tunneling. The situation becomes a bit complicated in the case of 
antiparallel alignment [6 = tt). When P = and 1, the spin current is zero, because for the 
former 1^ = ij_ while for the latter 7-j- = J| = 0, as shown in Fig. 4(b). For P ^ and 1, 
with increasing the bias voltage the spin current decreases steeply at the resonant position 
corresponding to the QD energy level Eq, and then keeps a constant until the other resonant 
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FIG. 5: The bias dependence of spin current I s in parallel (a) and antiparallel (b) configurations 
for different Coulomb interaction U, where P = 0.4, and the other parameters are taken the same 
as in Fig. 2. 

position corresponding to the level Eq + U where I s increases sharply, leading to the spin 
current shows a basin-like behavior. The larger the polarization P 0, 1), the deeper the 
bottom of the basin. These results demonstrate that the spin current I s has quite different 
behaviors in parallel and antiparallel configurations. 

Fig. 5 shows the bias dependence of spin current in collinear configurations for different 
Coulomb interaction U. It can be found that the overall qualitative behavior of I s against 
V looks similar to that shown in Fig. 4. With increasing U, in the parallel alignment the 
second plateau of the spin current J s becomes wider, as shown in Fig. 5(a); while in the 
antiparallel alignment apart from that the bottom of the "basin" becomes deeper, the width 
of the "basin" also becomes wider, as shown in Fig. 5(b). This is again from the fact that an 
increase of U makes the resonant positions corresponding QD discrete levels Eq and Eq + U 
become far separated. 
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FIG. 6: The current-induced spin transfer torque against the relative orientation angle for differ- 
ent polarization P (a) and different interaction U (b) at the bias voltage V = 25Fo/e. The other 
parameters are taken the same as in Fig. 2. 

C. Current-induced spin transfer torque 

The angular dependence of the CISTT for different polarization and Coulomb interaction 
is shown in Fig. 6. It is seen that the CISTT versus the angle 9 shows a sine-like behavior, 



consistent with the finding in Ref. [42|. At 9 = and tt, i.e., the magnetizations of the 
two ferromagnetic electrodes are collinearly aligned, the CISTT vanishes, which is obvious 
as the spin torque oc S# x (Sl x Sr), where Sl and Sr are the spin moments of the left 
and right ferromagnets. With increasing the polarization P, the magnitude of the CISTT 
is enhanced, as displayed in Fig. 6(a). This is in agreement with the statement that the 
torque is proportional to the polarization of the other ferromagnet 42 1. With the increase of 
Coulomb interaction U, the magnitude of the CISTT is decreased, as seen in Fig. 6(b). This 
can be understood in such a way that the contribution of the tunneling through the discrete 
level Eo + U is decreased with increasing the interaction U, thus leading to a decrease of the 
spin transfer torque. 

The bias dependences of the CISTT for different polarizations and Coulomb interactions 
at a given angle 9 = tt/3 are shown in Figs. 7(a) and (b). It is seen that with increasing the 
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FIG. 7: The bias dependence of the current-induced spin transfer torque for different polarization 
P (a) and U (b) at 9 = tt/3. The other parameters are taken the same as in Fig. 2. 

bias voltage, the CISTT first increases slowly to a sharp peak, then decreases dramatically 
to almost zero, and after undergoing a wider flat, it increases suddenly to another plateau. 
The position of the first sharp peak in the curve of CISTT versus V is observed to be 
independent of the polarization as well as the Coulomb interaction, showing that it is a 
resonant peak at the resonant position eq. The second sharp increase of the CISTT against 
V is also from the resonance at the QD discrete energy level Eq + U. From Fig. 7(a), 
one may find that with increasing the polarization the magnitudes of the CISTT at the 
two resonant positions increase. This is because the transmission coefficient of electrons is 

n 

proportional to the polarization [33J, while the CISTT depends on the tunneling electrical 
current that is determined by the transmission coefficient. Thus, an increase of polarization 
would enhance the CISTT at the resonant positions. The effect of the Coulomb interaction 
U on the CISTT is shown in Fig. 7(b). At the first peak, it is independent of U, suggesting 
that it is a resonance at the level Eo; while with increasing U, the second resonant positions 
move to higher voltages, showing that the second resonance takes place at the level Eq + U. 
It is interesting to note that at a given polarization the CISTT goes to a saturation plateau 
at the bias larger than (eq + U)/e, which is independent of the Coulomb interaction. From 
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FIG. 8: The current-induced spin transfer torque (a) and the spin current (b) against the electrical 
current as well as the current-induced spin transfer torque versus the spin current (c) for different 
U at = tt/3 and P = 0.4. The other parameters are taken the same as in Fig. 2. 

Fig. 7, we could see that the CISTT can be remarkably enhanced at the QD discrete energy 
levels which may be adjusted by changing the gate voltage. 

What is the relationship between the CISTT, the electrical current and the spin current 
in this FM-QD-FM coupled system? The answer is presented in Fig. 8 for different U at 
9 = tt/3. With increasing the electrical current, the CISTT goes up linearly to a round 
peak, then decreases slowly to a broad minimum, and then increases sharply, as shown in 
Fig. 8(a). The first round peak is owing to the resonant tunneling at the level £q] while 
the second sharp increase is caused by the resonant tunneling at the level Eq + U. The 
spin current I s is found to increase non-monotonically with increasing the electrical current 
/, as illustrated in Fig. 8(b). This is quite different from the spin current in FM-SC-FM 
tunnel junctions where the spin current is proportional to the injection current [43]. In the 
curves of I s against /, the two resonances can be clearly seen. The CISTT as a function 
of I s shows a kink-like behavior, as displayed in Fig. 8(c). With increasing I s , the CISTT 
first shows a linear behavior, then suddenly drops to nearly zero, and then increases almost 
linearly again. This kink-like behavior is caused by the resonant tunneling at the levels Eq 
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and Eq + U, which manifests the same thing as in Fig. 7. Away from the kink, the CISTT 
appears to be proportional to the spin current. The reason is that the electrons with opposite 
spins from the left ferromagnetic electrode entering into the right electrode exert the torques 
with opposite directions on the spins of the right electrode, giving rise to that the magnitude 
of the CISTT is determined by the difference between the magnitudes of torques exerted by 
electrons with spin up and spin down. Considering that the spin current is defined by the 
difference between the electrical currents of spin up and spin down, while the spin transfer 
torque is proportional to the electrical current jf], Q], the proportionality between the CISTT 
and the spin current is conceivable. It is seen that for different U the curves in Fig. 8 are 
almost unchanged, showing that the Coulomb interaction has less effect on the CISTT versus 
/ and I s , as well as I s versus /. This is because the Coulomb interaction U only changes 
the resonant positions. 

IV. SUMMARY 

We have investigated the spin- dependent transport properties in a FM-QD-FM coupled 
system by means of the nonequilibrium Green functions. It has been found that by changing 
the relative orientation of both magnetizations, the spin current can change its magnitude 
even its sign. For positive bias voltages, when the relative orientation angle is less than 7r/2, 
the spin current shows a step-like behavior, and the spin current is positive; when the relative 
orientation angle is greater than n/2, the spin current behaviors similar to a basin-like shape, 
and is negative at the bottom of the "basin". With increasing the Coulomb interaction U in 
the QD, it has been uncovered that the resonant positions for the tunneling electrical current 
and spin current become far separated. The CISTT is observed to first increase slowly to 
a sharp peak with increasing the bias voltage, then decrease dramatically to almost zero; 
and after undergoing an unchanged stage, it increases suddenly to another plateau. Such 
a behavior is obviously resulted from the resonant tunneling through the central QD. The 
magnitudes of the CISTT at two resonant positions are observed to increase with increasing 
the polarization of the ferromagnet. At the bias larger than (e + U)/e, the CISTT is seen 
to reach a saturation plateau which is independent of the Coulomb interaction. In addition, 
it has been demonstrated that the CISTT as a function of the spin current shows a kink-like 
behavior, and away from the kink, the CISTT increases almost linearly with increasing the 
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spin current. The CISTT as a function of the electrical current shows a non-monotonic 
behavior, while the spin current is found to be nonlinearly proportional to the electrical 
current. Besides, the Coulomb interaction is shown to have less effect on the behaviors of 
the CISTT as functions of the tunneling electrical current and spin current, as well as the 
spin current as a function of the electrical current. 

Finally, we would like to mention that the CISTT can be greatly enhanced when the 
bias voltage meets with the discrete levels of the QD at two resonant positions. While the 
energy level of the QD can be adjusted by tuning the gate voltage, the resonant property of 
the CISTT offers an alternative premise to develop a spintronic device through the current- 
controlled magnetization reversal effect. 
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